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We study the formation of central naked singularities in spherical dust collapse with a cosmological 
constant. We find that the central curvature singularity is locally naked, Tipler strong, and generic, 
in the sense that it forms from a non-zero-measure set of regular initial data. We also find that the 
Weyl and Ricci curvature scalars diverge at the singularity, with the former dominating over the 
latter, thereby signaling the non-local origin of the singularity. 
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I. INTRODUCTION 

One of the outstanding issues in general relativ- 
ity is that of the final state of gravitational collapse, 
a crucial aspect of which is the possibility of space- 
time singularities — events at which Riemannian curva- 
ture typically diverges, the spacetime is geodesically in- 
complete, and any classical theory of gravity necessarily 
breaks down. It has long been known that under a va- 
riety of circumstances, spacetimes which are solutions of 
Einstein's equations with physically reasonable regular 
initial data, inevitably develop singularities 

In an effort to protect the applicability of general rel- 
ativity, Penrose conjectured that such singularities must 
be hidden by an event horizon, and thus invisible to an 
asymptotic observer, i.e., they cannot be globally naked 
[^|. This constitutes in essence what has become known 
as the weak cosmic censorship conjecture. However, it is 
quite possible — at least in principle — for an observer to 
penetrate the event horizon and live a rather normal life 
inside a black hole. This motivated the strong cosmic 
censorship conjecture, which broadly states that timclike 
singularities cannot occur in nature, i.e., there are no 
locally naked singularities || . 

A lack of tools to handle global properties of the Ein- 
stein equations (and respective solutions), together with 
their high non-linearity, have been the main obstacle to 
provable formulations of the (weak or strong) cosmic cen- 
sorship conjecture. Whilst efforts are being undertaken 
in this direction Q , one can hope that the detailed study 
of specific models helps to isolate some defining features 
of singularity formation and structure, thereby contribut- 
ing towards a precise, counter-example-free formulation 
of the conjecture. 

One such model is inhomogeneous Tolman-Bondi col- 
lapse HU , whose general solution is analytically obtain- 
able in closed (or parametric) form. For this reason, 
and because spherical symmetry is, arguably, a plausi- 
ble approximation for the geometry at the late stages 
of collapse 0, Tolman-Bondi collapse has been exten- 
sively studied by many authors The main results 
stemming from these analyses, pertinent to the context 



of singularities and cosmic censorship, may be summa- 
rized as follows: (i) A central curvature singularity al- 
ways forms; (ii) This singularity can be either locally or 
globally naked, depending on the initial data; (iii) The 
singularity is generic, in the sense that an infinite number 
of outgoing non-spacelike geodesies terminates at the sin- 
gularity in the past; (iv) The singularity is gravitationally 
strong in the sense of Tipler . 

We remark that inhomogeneity plays a key role in the 
global visibility of the singularity in asymptotically flat 
spherical dust collapse, with the former being uniquely 
determined by the derivatives of the initial central energy 
density profile. In connection with the cosmic censorship 
conjecture and the physical meaningfulness of the singu- 
larity, it should be noted that the central naked singu- 
larity of inhomogeneous Tolman-Bondi collapse is always 
Tipler strong, irrespective of the initial data and, 
furthermore, it is marginally stable against non-spherical 
linear perturbations Jl2[ . 

Recent observations of high-redshift Type la super- 
novae |[3|,|l4| and peculiar motion of low-redshift galaxies 
]l5| , appear to indicate that the present radius of the uni- 
verse is accelerating, thus suggesting the existence of a 
positive cosmological constant, A > 0. This has sparked 
a renewed interest in gravitational collapse with a cos- 
mological constant (TfJ. Markovic and Shapiro |l7| ana- 
lyzed the spherical homogeneous collapse of a dust cloud 
in the presence of a positive cosmological constant, and 
found that depending on the initial data (cosmological 
constant, gravitational mass of the cloud, and its comov- 
ing radius), the resulting spacetime can be either (i) a 
Schwarzschild-de Sitter black hole |lS[| , (ii) a bouncing 
sphere, or (iii) a de Sitter-like global cosmological singu- 
larity. Their analysis was qualitatively generalized to the 
inhomogeneous and degenerate cases, for both A > and 
A < 0, by Lake The collapse of null dust with a neg- 
ative cosmological constant was studied by Lemos [ pfj[ , 
who showed that a Tipler strong, globally naked singu- 
larity develops for spherical collapse (but not for toroidal, 
cylindrical, or planar geometries). 

In this paper, we examine in detail the effects of a pos- 
itive cosmological constant in the singularity formation 
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and structure in spherical dust collapse in an asymptoti- 
cally de Sitter spacetime. The existence of A > changes 
the nature of the solution not only at large radii, but also 
near the center, and thus the central singularity if it 
exists — may have a different structure from that occur- 
ring in asymptotically flat spherical dust collapse. The 
current status of the cosmological constant makes perti- 
nent the analysis of its effects on singularity formation 
and structure, particularly in the cosmic censorship con- 
text. Does a central curvature singularity exist in spher- 
ical dust collapse with A > 0? If it does, is it naked and 
gravitationally strong? As we shall show below, both 
questions are answered positively. 

For definiteness, we shall hereafter refer to spherical 
dust solutions of the Einstein equations with a cosmo- 
logical constant as Tolman-Bondi-de Sitter spacetimes. 

The paper is organized as follows: Section II derives 
the Tolman-Bondi-de Sitter metrics from the Einstein 
equations. In Sec. Ill, the existence and genericity of 
a central curvature singularity are proven. Section IV 
discusses the singularity's curvature strength and the rel- 
ative strength of the Weyl and Ricci scalars. Section V 
concludes with a summary and discussion. 

Geometrized units, in which G = c = 1, are used 
throughout. 



Introducing the auxiliary functions 
k(t,r) 



1 - e 2 *i?' 2 , 



m(t, r) 



1 



-R R 



f?2 



Einstein's equations^ simplify greatly to 

R 2 = 2mRT 1 -k + -R 2 , 
3 

k = 0, 
m = 0, 



(8) 
(9) 

(10) 

(11) 
(12) 

(13) 



with the constraint 

m = 4irR 2 R'p(t,r) 
The Tolman-Bondi-de Sitter metrics are then given by 
R' 2 



ds z = -dt 2 



1 - k 



dr 2 +R 2 dtt 2 , 



(14) 



where R(t : r) is a solution of Eq. (|T(j), with initial data 
given by Eq. (pi), 



m(r) = 4tt / i? 2 (0, r)R'(0, f)p(0, f)df, 



(15) 



II. TOLMAN-BONDI-DE SITTER SPACETIMES 

The Tolman-Bondi-de Sitter family of solutions is 
given by a spherically symmetric metric, written here in 
normal Gaussian coordinates {t, r, 9, </>}: 



ds 2 = -dt 2 + e-^^dr 2 
dn = d0 2 + sin 2 P^ 2 



R 2 (t, r)dtt 2 



together with the stress-energy tensor: 

Tab = p(t, r)u a u b - -^-g a b = pS a Sl - -^-g a b, 

07T 07T 



(1) 

(2) 



(3) 



where u a — Sf is the 4-velocity of a dust element, p(t,r) 
is the energy density, and A the cosmological constant. 

With the metric (|l|), the independent non- vanishing 
Einstein tensor components are 



G u = R- 2 [-Re 2 *(2R'm' + 2R" + R~ R ) 



G r t 

Gag 



-2RVR+1+R 2 }, 
-2R- l (R' + R'm, 



= —R~ 



sin 



e --(2i?i? + i? 2 + l)-.R' 2 
G^ = R{M> + i?'*'e 2 * 



+i?"e 2 * -R+$R- 4/ 2 R), 



(4) 
(5) 

(6) 
(7) 



where the overdot and prime denote partial differentia- 
tion with respect to t and r, respectively. 



together with an initial velocity profile R(0, r) [which 
fixes k(r), via Eq.(^FJ) evaluated at t = 0], 

Without loss of generality, we consider here the k = 
case (corresponding to gravitationally unbound matter 
configurations), since it allows for an analytical solution 
of Eq. (10) in closed form: 



R&r) =(^)"sinh 2 / 3 T(t,r), 



T(t,r) 



'3A 



[*c(r) - 1] , 



(16) 
(17) 



where t c (r) is the proper time for complete collapse 
[R(t c ,r) — 0] of a shell with initial area radius R(0,r), 
which is fixed by Eq. (n6J) at t = 0: 



(18) 



3A 



sinh 




where the scaling R(Q,r) = r was adopted [note that, 
when k = 0, R(0,r) is automatically fixed by the choice 



* Since there are only three functions to be determined and 
four equations, only three of these are independent, with the 
remaining one acting as a constraint. We take Eqs. (|Io[)-(p^) 
as our complete set, and Eq. ( |l3| ) as the constraint equation, 
since it provides a simple relation between the initial data and 
the initial mass profile. 
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for the radial gauge via Eq. (|10|)]. The relevant deriva- 
tives of the area radius are 



R(t,r) = -y^i?cothT, 
where the minus sign corresponds to implosion. 



(19) 
(20) 



III. 



EXISTENCE AND GENERICITY OF 
SINGULARITIES 



From Eq. @, we have R = ± v / 2mR- 1 + {A/3)R 2 , 
where the plus or minus sign corresponds to expansion or 
implosion, respectively. If R < 0, then every dust shell 
implodes and inevitably collapses to vanishing proper 
area in a proper time given by Eq. (|l8|). If R > 0, pro- 
vided the initial "acceleration", i?(0, r) = —^jfi- + jr < 
0, all the shells will initially expand towards increasing 
area radius, reach a maximum value R max (r), and then 
collapse back through their original radii, eventually end- 
ing up with zero proper area. In either case we have 
R(t c (r),r)=0. 

It then follows that at t = t c (r) the Kretschmann 
scalar, 



xr — d Tyabcd o 

A- = Jriabcdtt — o 



ml mm' m' 2 

'r^~ 8 r^ + 12 1v> (21) 



diverges, thereby signaling the existence of a curvature 
singularity. It has been shown (explicitly for the A = 
case, but the result also holds for A > 0) that the cur- 
vature singularity for events with r > is spacelike, and 
thus physically uninteresting . Of potential interest is 
the central curvature singularity that forms at the event 
(t = t c (0),r = 0). 

We want to determine if such a singularity is (at least 
locally) naked, i.e., if there exists at least one future- 
oriented radial null geodesic with past endpoint at the 
singularity. To do so, we examine the outgoing radial 
null geodesies equation: 

where Eqs. (0) and (||) were used. Expanding p(0, r) = 
p c {f) near r = 0, 



i=0 



(23) 



near the singularity (where r, T — > + ) we have, to lead- 
ing order, 



m{r) = m r 3 + m n r n+3 + 0(r n+4 ), 
*c(r) =t + t n r n + O(r n+1 ), 

2) r° 

+ 0{r n+2 ) x 0(Ti), 



(24) 
(25) 



R(t,r 



m$r + M n r n+1 )(t + t n r n -t)z 



(26) 



where t n and M n are real coefficients linear in m n = 
(4ir/n)p n , with n > 0; p„ = (d n p c /dr n ) r =o is the first 
non- vanishing derivative of the central energy density dis- 
tribution, and 



to = sinh 1 



/ A 
6mg 



(27) 



We now follow the method outlined by Barve et al. 
p2|, and assume that there is a regular solution to the 
outgoing null geodesies equation near r = 0. To leading 
order in r, we write such a solution as 



t = t + ar a , 



(28) 



where a, a £ M + . We note that, since T > 0, we must re- 
quire a > 71. If a = n, we have the additional constraint, 
a < t n . 

From Eqs. © and ©, we obtain 



/7 ( ,. ri : ('^£yi| r ¥+i +0 (r-+ 2 - 



(29) 



Let us first consider the a > n case. Differentiating 
Eqs. (H) and (f2§) w.r.t. r, we obtain, from Eq. (||), to 
leading order, 



, , < 2n 

acrr*- 1 = f y + 1 



Self-consistency fixes 



9m n \ 3 1 



2n 

/9m o y s 
«=( — ) *i 



(30) 

(31) 
(32) 



The condition a > n now reads n < 3. For n = 1,2 
(i.e., for pi ^ 0, or pi = and p2 ^ 0) there is a self- 
consistent solution to the outgoing radial null geodesies 
equation in the limit t — > to, r — > 0, and thus there is at 
least one outgoing radial null geodesic starting from the 
singularity, which is therefore naked. 

We now examine the case n = a = 3. Proceeding as 
before, we obtain, to leading order, 



3ar z = 3 



9mr 



(£ 3 -a) 3 r 2 , 



(33) 



which is identically satisfied provided 
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M 3 a 2 + 2M%a = 0, 



(34) 



where M = (9mo/2) 1 ' 3 . This equation has two non-zero 
distinct roots (other than the a = trivial root), given 

by a = (M 2 /2) ± VM 4 - 8i 3 , if t 3 < ± (^f 1 )', which 
imposes a constraint on /?3, for a given po- ln addition, 
self-consistency also requires that a < t%, which leads to 



^-4 +V ^ 



1 



16 - 4M 2 <U< -M 



(35) 



Thus, as long as one restricts ourselves to initial data that 
satisfies the above condition, the singularity is naked. We 
note, however, that the n — 3 case is a less generic case 
than the n < 3 one, as it requires that p\ = p2 = and 
ps obey condition (|35|). 

Let us now investigate whether there is only one null 
geodesic emanating from the singularity (in which case 
the singularity would be "visible" for an infinitesimal 
amount of time), or an entire family — in which case the 
singularity would be visible for an infinite amount of 
time. We write the equation for the outgoing radial null 
geodesies to next order as 



t = t + ar a + br° 



(36) 



where a, a, S 6 R + . Again, we shall consider first the case 
a < n. Proceeding as before, we have 

aar"- 1 + b{a + < 5)r CT+<5 - 1 = Ar^ a + Br 5+ ^, (37) 

where 



, s (f +1 )(^y ( j, 



(38) 
(39) 

(40) 

thus fixing 6 = 1 - f , and b = -§ (Sffla)' t~*. Hence, 
to this order there is a single future-oriented radial null 
geodesic emanating from the singularity. 

Let us now examine the case a = n = 3. From Eqs. 
(EQ) and (pq), we get 



Self-consistency fixes a and a as before, and implies 
6^ + 1 + ^ r 5-i = Br -f ) 



3+<5 



where 



R{t,r) = air 3 + a 2 r 



/9m \3 5 
"!=( — ) (*a-o)< 



"2 = -gO:i(*3 - a)' 



(41) 

(42) 
(43) 



Taking the partial derivative of R(t,r) w.r.t. r and 
equating it to that of Eq. (Bf?) yields 



3ar 2 + 6(3 + 6)r 



2+8 



3air 2 + a 2 (3 + 5)br 



2+5 



(44) 



Self-consistency at first-order fixes a = a\ [which is 
equivalent to the cubic equation (Q)] as before, and a 
solution exists if condition (|35| ) is satisfied. At the next 
order, b drops out and self-consistency requires a.% = 1. 
This amounts to a particular value of a, which is consis- 
tent with those determined from Eq. (|3^), if and only 
if 



ta 



M 2 i 

M H — 4 ± VM 4 - 8M 3 + 4M 2 + 4. (45) 



Positivity of the radicand requires M < 1.076381196, or 
M > 7.453694050, but the first range is disallowed since 
we must have £3 > 0. Consistency between condition 
© and Eq. © further imposes M > 8.799799016. 
Hence, a fully self-consistent solution with n — 3 exists 
when M > 8.799799016 and t 3 is given by Eq. @. 
In such case, an entire one-parameter family of outgoing 
radial null geodesies (parameterized by b) departs from 
the singularity. Clearly, this is a very special case — it is 
a set of measure zero in the initial data|] — and in general 
there is only a single outgoing null geodesic with past 
endpoint at the singularity. 



IV. GLOBAL VISIBILITY 

We have shown that there is a generic central curva- 
ture singularity in Tolman-Bondi-de Sitter collapse. This 
singularity is at least locally naked. It would be globally 
naked if outgoing geodesies starting from the central sin- 
gularity at (t — to, r — 0) could reach future null infinity. 
In the present case, this would require future-directed 
geodesies to avoid the event horizon (i.e., remain outside 
it at all times) and cross the cosmological horizon. 

Let us first examine the possibility of event horizon 
avoidance. In spherical dust collapse, the event hori- 
zon (EH) coincides with the apparent horizon (AH) at 
the boundary of the spherical mass distribution, r = r m . 
Hence, in order for the outgoing null geodesies to escape 
the EH, they must cross r = r m before the AH, thereby 
avoiding becoming trapped, and hence eventually ingo- 
ing. 

The AH is a spacelike 2-surface defined by the locus of 
events where null wavefronts become "frozen" , and in the 
adopted coordinates it is given by R :a R,bg ah = 0. With 
the metric (|l4| ) and Eq. ([To]) we have then 



We note that, while this is obvious in the present case, 
it is in general unclear what measure or topology should be 
imposed on the space of initial data. 
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2m = 0. 



(46) 



This equation has three distinct real roots if 3m\/A < 1 , 
two of which are positive and given by 



i?i = — = sin 
VA 



R-2 



sin 



with i?2 > 
-i 



2 

Ri > 



1 

-sm 



-sm 



L (3m\/A) 



L (3m7A) + — 



(47) 



(48) 



0, corresponding to the choice < 
sin -/3 < tt/2, < /3 < 1. The third root, R 3 = 
—R1 — R2, is negative and hence unphysical. R2 is a gen- 
eralized cosmological horizon (R2 = when m = 0) 
and i?i the black hole apparent horizon (i?i = 2m when 
A = 0; the apparent and event horizons coincide in the 
static case). For 3m\/A = 1, the two horizons coincide. 
If 3r7i\/A > 1, there is one negative real root and two 
complex (conjugate) roots, all of which are unphysical: 
the spacetime does not admit any horizons in this case. 
FromEq. @ together with Eqs. we obtain 



tAu{r) = t c (r) - -j= sinh 1 



p \ a 



(p sin 9) 2 



where 



a 



e 



2 

6m 

p 

1 . 

— sin 
3 



. (49) 



(50) 
(51) 



(52) 



A necessary condition for the singularity to be glob- 
ally naked is ^ah(O) > t c (0), which ensures that ra- 
dial null geodesies emanating from the singularity do 
so before the AH forms. However, from Eq. ( [49] ) 
we have ^ah (0) = t c (0) - X, where < X < 
(2/V3A)sinh~ 1 [(V3/m) 1 / 2 /A]. Outgoing radial null 
geodesies are therefore unavoidably trapped inside the 
AH and the singularity cannot be globally naked. 



Dust cloud 
surf Lice 



Singularity 




FIG. 1. Central naked singularity in spherical Tol- 
man-Bondi-de Sitter collapse. The singularity is locally 
naked: it forms before the boundary of the mass distribu- 
tion undergoes complete collapse, but after the apparent (and 
thus event) horizon form. Null geodesies emanating from the 
singularity are necessarily trapped. 



V. CURVATURE STRENGTH 

A crucial property of a singularity is its curvature 
strength. A singularity is said to be gravitationally 
strong in the sense of Tipler |l0[ if every collapsing vol- 
ume element is crushed to zero at the singularity, and 
weak otherwise (i.e., if it remains finite). It is generally 
believed — although not yet proven p3[ — that spacetime 
is geodesically incomplete at a strong singularity, but ex- 
tendible through a weak one |l(],^4|,^5| . 

A precise characterization of Tipler strong singularities 
has been given by Clarke and Krolak (26) , who proposed 
(among other conditions) the strong focusing condition: 
There is at least one null geodesic, with tangent k a and 
affine parameter A (with A = at the singularity), along 
which the following is satisfied: 



lim X 2 R ab k a k" > 0. 



(53) 



This is a sufficient condition for the singularity to be 
Tipler strong and corresponds to the vanishing of any 
two- form defined along such a geodesic, at the singularity, 
due to unbounded curvature growth. 

Let us now consider a radial null geodesic with tangent 
k a = (fc*,fc r ,0,0), where fc° = and 



fc* = FR = R'k r , 



(54) 



where F can be written as an explicit function of the 
affine parameter, F = F(X), obeying the differential 
equation (which follows from the geodesic equation, 
k a V a k b = 0): 



dF 



F 2 1 + R + R 



R!_ 
R' 



(55) 



From Eqs. (|l4|) and (i^0|) we have then 



2(k 



t\2 



R?R r 
(56) 



Using Eqs. and @-(||), we get 



lim X 2 n = lim X 2 F 2 — 

A^O A^O R 



m 2 m 
— = urn y — — 

t^t ,r^o A R' 



K n lim x 



(57) 
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where l'Hopital's rule was used twice in the second equal- 
ity, and 



X{t,r) 




4m A p2 _ o rti' 

Y _ R 3 11 Z R' 



where Eq. ( |io| ) was used, and 

m! 1 f 2n \ ,. m 
K n = hm — = hi hm — 



(58) 
(59) 



K lim r 2 = 0, 



2m^\ s 



for n < 3, (60) 



#3 = [^j (*3 - a) 5 > 0, torn - 3, (61) 

X = 3fy + lJ (6m 2 ) 3 t„ 3 , forn<3. (62) 

[Note that K% exists as long as the initial data obeys 
condition ([35])]. 

Therefore, for n < 3, 



lim y = lim 1 - f(ri)y Kr 



3— n 



1, (63) 



where /(rc) = [(8n/3) - (2/3)]/ x /2[(2n/3) + l]/3 > 0. 
If n = 3, we have 



lim x 2 - (1 - V2^V 



t->to,i"->0 



- 2 = v 2 

— Xo> 



(64) 



which is positive definite provided K3 ^ 1/2. 

Summarizing, the singularity is gravitationally strong 
in the sense of Tipler if n = 3. For n < 3, 
lirriA— >o X 2 R a bk a k b — 0, and the singularity may or may 
not be Tipler strong. 



A. Weyl curvature 

In order to gain further insight on the nature of the 
central singularity, we examine the behavior of the Weyl 
curvature scalar in the limit t — > to at r = 0. In a four- 
dimensional manifold, the Weyl tensor is defined as 

Cabcd = Rabcd + g a [dRc]b + 9b[ C Rd]a + g R-9a [c9d] b , (65) 

and represents the part of Riemannian curvature which 
is not locally determined by the matter distribution. A 
divergence of the Weyl scalar, C = C a bcdC abcd , indicates 
a blow-up in curvature caused not by the local matter 
distribution — as in the case of the Ricci scalar, or any 
other scalar constructed solely from the Ricci tensor — 
but by the matter content of the spacetime at other 
points. 



In the metric (JIJ), the non- vanishing independent com- 
ponents of the Weyl tensor are 



Cfrtr — 2R'' 



Ctete — 



m \ 



R 2 R> R 3 J 



R 3R'' 



C r 8r8 — 



m'R' mR' 



3 

• 2 , 



R 



Cr4>r<f> — Sin 0C r g r 

Com = 2 sin 2 ( mR - 



m'R' 
Wf 



The Weyl scalar is 
C(t,r) 



R 4 



m 
3R 1 



rn 
~R 



(66) 

(67) 
(68) 
(69) 
(70) 
(71) 



(72) 



At the singularity we obtain, for n < 3, 

C si ng = lim C(t,r) = C lim r^ 1 ^ = +oo, (73) 

t— >to,r— >0 r^0 

where C* = 48(2Xn/9) 2 (9m /2)- 4 / 3 i I 7 8/3 > 0. 
If n = 3, 

- hm 48 (^) t^Klr- u = +oo. (74) 

Thus, for all the initial data leading to naked singulari- 
ties (n < 3), the Weyl scalar diverges along outgoing null 
geodesies at the singularity. We note that the present 
analysis, for outgoing radial null geodesies, remains un- 
changed up to a sign for dt/dr and a in Eq. (^2|), for 
ingoing radial null geodesies. We conclude, therefore, 
that the Weyl curvature scalar diverges at the singular- 
ity along both outgoing and ingoing null geodesies. This 
is in agreement with the results of Barve and Singh pjj ] 
for asymptotically flat spherical dust collapse, and goes 
against the speculation by Penrose ]23] that the Weyl cur- 
vature should diverge along ingoing geodesies, and vanish 
along outgoing geodesies terminating at the singularity. 

Finally, we compute the Ricci curvature scalar, 1Z, at 
the singularity and compare it to the Weyl scalar, to 
determine the relative contributions of curvature (i.e., 
locally versus non- locally induced). The Ricci scalar is 



2m' 
R 2 R' 



4A. 



At the singularity we have, for n < 3, 
ft sing = lim Tl(t, r) = TZ Q lim r^ 2 " 



where TZ = 2(9m /2)- 2 / 3 t„ 4/3 > 0. 
For n = 3 



(75) 



(76) 



6 



r— *0 \ 



2 / 



" 3 r 6 = +00. 



(77) 



Hence, {C /1Z) s \ ng diverges as r - 10 ™/ 3 ; for n < 3, and as 
r~ 6 for rt = 3. Such a dominance of the Weyl curvature 
over Ricci indicates a predominantly non-local origin for 
the Riemannian curvature divergence at the singularity. 



VI. CONCLUSIONS 

We have examined in detail the central curvature sin- 
gularity occurring in general inhomogeneous spherical 
dust collapse in an asymptotically de Sitter spacetime. 
This singularity was found to be locally naked for a wide 
class of initial data: p\ ^ 0, or p\ =0 and pi ^ 0, or 
pi = P2 = and P3 ^ 0. In the latter case (if, in ad- 
dition, P3 is determined by po), an entire one-parameter 
family of outgoing radial null geodesies has its past end- 
point at the singularity, which is thence visible for a finite 
amount of time. One must note, however, that this case 
is non-generic, whereas the first two are generic insofar 
as initial data is concerned. Hence, generically, a single 
null geodesic escapes the singularity. 

Irrespective of the initial data, the singularity cannot 
be globally naked, since it forms after the apparent hori- 
zon does, and, consequently, any null geodesies emanat- 
ing from it are necessarily trapped and thus cannot es- 
cape the event horizon. 

Regarding curvature strength, our analysis, based on 
the sufficient condition by Clarke and Krolak po| , re- 
vealed that the singularity is Tipler strong for the case 
n = 3 and may or may not be Tipler strong if n < 3. We 
note that, while more restrictive than the cases n < 3, 
the n = 3 case is of finite measure in the space of initial 
data — P3 has a finite real range, given implicitly by con- 
dition (|35| ) — and hence generic, in the sense that it has 
codimension two in the (countable, infinite-dimensional) 
space of derivatives of p c {r). 

We also found that both Weyl and Ricci scalars diverge 
at the singularity along ingoing and outgoing geodesies. 
The Weyl scalar divergence dominates over the Ricci, in- 
dicating a predominantly non-local origin of the Rieman- 
nian curvature unbounded growth. Interestingly, while 
features such as global visibility and curvature strength 
appear to depend critically on the local matter distri- 
bution, the divergence of Riemann curvature invariants 
(e.g., Kretschmann scalar) seems to be associated with 
the matter distribution at other points. It is tempting to 
speculate that one could alter the local matter distribu- 
tion, so as to cure such "problems" as local visibility and 
Tipler strong strength, while still maintaining a diver- 
gent Riemannian curvature invariant (built solely from 
the Riemann tensor, without any of its internal index 
contractions). This naturally leads to the question of 
whether spacetime can be extendible through a Tipler 



weak singularity where curvature is, nevertheless, diver- 
gent. 

Finally, in the context of cosmic censorship, with the 
current status of A, any realistic formulation (of the 
strong version, at least) might have to exclude the cases 
corresponding to A > 0, if a stability analysis — to be 
defined in a suitable, precise manner — would reveal the 
singularity to be a persistent (i.e., exist for all times, with 
the same properties) feature of the spacetime. The re- 
sults of Deshingkar, Joshi, and Dwivedi [jnj, and those 
of Harada, Iguchi, and Nakao [|l2| , reveal that the central 
singularity in spherical dust collapse with A = is sta- 
ble (i.e., remains locally naked and Tipler strong) against 
initial data perturbations, and marginally stable against 
metric (and matter coupled to metric) perturbations. The 
similarity between such a singularity and the one dis- 
cussed in the present paper — the conditions for visibility 
and the slope of the outgoing radial null geodesies [cf. 
Eq. (|3l|)1 with past endpoint at the singularity are the 
same in both cases (see, e.g., p3]) — suggests that anal- 
ogous stability properties may also hold for the central 
curvature singularity in Tolman-Bondi-de Sitter collapse. 
This issue is currently under investigation p9| . 
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